Momentum transfer to small particles by aloof electron beams 
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The force exerted on nanoparticles and atomic clusters by fast passing electrons like those em- 
ployed in transmission electron microscopes are calculated and integrated over time to yield the 
momentum transferred from the electrons to the particles. Numerical results are offered for metallic 
and dielectric particles of different sizes (0-500 nm in diameter) as well as for carbon nanoclusters. 
Results for both linear and angular momentum transfers are presented. For the electron beam 
currents commonly employed in electron microscopes, the time-averaged forces are shown to be 
comparable in magnitude to laser-induced forces in optical tweezers. This opens up the possibility 
to study optically-trapped particles inside transmission electron microscopes. 
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I. INTRODUCTION 

Electromagnetic forces in optical tweezers are currently 
employed to trap small particles ranging in size from 
nanometers to several microns 0, Q, and to manipu- 
late them in all spatial directions This type of 
forces is also used to characterize the elastic properties 
of deformable tiny objects (e.g., living cells (3), to ob- 
tain quantitative information on mechanical properties 
at small length scales Q, and in general, to fix the posi- 
tion of those particles so that they can be manipulated 
at will. 

In this context, transmission electron microscopy offers 
a potentially useful tool to study optically trapped par- 
ticles, providing excellent spatial resolution (sometimes 
below 1 A) when sub-nanometer electron beams are em- 
ployed [g, while allowing spectroscopic characterization 
with sub-eV accuracy. Actually, transmission electron 
microscopes are routinely exploited to probe local op- 
tical response properties @, and more recently, also to 
determine photonic structures of complex materials Q . 

A major problem that may arise when combining elec- 
tron microscopy with optical tweezers or other types of 
optical trapping (e.g., optical lattices [HEIE]]) is that 
the passing electrons can kick the particles out of the 
trapping locations (see Fig. QJ. In this work, we show 
that the momentum transferred from the passing elec- 
trons to the particles can be well below the threshold 
needed to kick them out for commonly employed trapping 
laser intensities, although a detailed comparison between 
trapping forces and electron-induced forces suggests that 
both weak and strong perturbation regimes are possible 
depending on the distance between the particles and the 
beam, all of them within the range that allows a suf- 
ficiently large electron-particle interaction as to perform 
electron energy loss spectroscopy (EELS) with significant 
statistics for in vacuo optically-trapped particles. 

The moving electrons can be in fact regarded as a 
source of evanescent electromagnetic field that probes the 
sample locally, and in this sense, they can be also used 
to produce deformation in elastic particles, oscillations 
of trapped particles around their equilibrium positions, 
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FIG. 1: (color online). Schematic representation of the pro- 
cess considered in this work: a fast electron moving with im- 
pact parameter b and velocity v with respect to a polarizable 
particle transfers momentum Ap = (Ap x , Ap z ) to the particle 
via electromagnetic interaction. 



and other interesting effects associated to the transfer of 
momentum within accurately controlled spatial regions. 



II. THEORY 

The electromagnetic force exerted on a particle in vac- 
uum is given by the integral of Maxwell's stress tensor 
over a surface S embedding the particle 01 as 



F(i) = — / ds[E(s,t)E(s,i) • n + H(s, i) H(s, i) • n 
|E(s,t)| 2 + |H(s,i)| 2 )], 
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where n is the surface normal and Gaussian units are 
used. The momentum transferred to the particle, Ap, is 
obtained by integrating of F(t) over the time. This yields 



Ap = jF{t)dt = J^ F{cj)duj, 



(1) 
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where 



F(w) = 4^M/ dB[E( B) w)(E(s,w)-n)* 
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+H(s,w) (H(s,w)* • n)* 
-|(|E( 8>W )| 3 + |H(8 >W )| a )]}, 

and the Fourier transform is defined as E(r, w) = 
/ dtE(r,t)exp{iwt}. 

The force acting on the particle is due in part to radi- 
ation emitted as a result of interaction with the electron 
and in part to the reaction force experienced by the pro- 
jectile. For small particles, the effect of radiation emis- 
sion is negligible and the trajectory is deflected by an 
angle f=a Ap/mv, where m and v are the mass and veloc- 
ity of the electron. Non-retarded calculations have shown 
that this angle is too small to be easily measured fl3|. 



A. Small particles 

Let us first consider a small isotropic particle suffi- 
ciently far away from the electron beam as to neglect 
higher multipoles beyond induced dipoles. The particle 
is then characterized by its frequency-dependent polar- 
izability cx(u>), and the force exerted by each frequency 
component of the external field E(r, uS) reduces to [Ty 

F(w) = Re{a ^ J Bf t (r,w)V[^ cxt (i-,w)]*}- (3) 



This expression can be derived from Eq. J3J) by consider- 
ing an integration surface arbitrarily close to the object 
and by using the expressions for the electric and magnetic 
fields induced by a small polarizable particle in terms of 
its polarizability a. For an electron moving with velocity 
v towards to positive z direction and passing by the ori- 
gin at t = 0, the external field is readily calculated from 
Maxwell's equations to yield 
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where R = (x,y) and 7 = 1/ yl — v 2 / c 2 . Inserting Eq. 
into Eq. one obtains 



Fu > = %^[-Re{a} /'(— ) x + 2 7 lm{a} /(— ) z],(5) 



where 



f(C) = K 2 (C)+K 2 (()/ 1 2 



(-6,0) 



and the particle is taken to be situated at R 
with respect to the beam (see Fig. Q). 

Symmetry considerations lead to the conclusion that 
Rayleigh scattering of the external-electron evanescent 



field Q| produces a radiation pattern with inversion sym- 
metry with respect to a plane perpendicular to the tra- 
jectory. This means that the overall transfer of momen- 
tum to the induced radiation is zero in the small-particle 
limit, so that Ap z accounts for all momentum transfer 
to the moving electron along z. Then, the contribution 
of each lu component to the electron energy loss rate is, 
within the non-recoil approximation valid for sufficiently 
energetic electrons, vF z (uj). Actually, one finds that the 
identity vF z (uj) = HujP(uj) is satisfied, where P(w) is 
the frequency-resolved loss probability as previously ob- 
tained for small particles As a consequence, F z van- 
ishes in the u) — * limit, since P(uj) remains finite. 

This behavior is quite different from F x , which goes 
to a finite value for small oj's, namely F x {oj = 0) = 
4e 2 Re{a(0)}/u 2 6 3 . (Incidentally, momentum transfer 
along x produces negligible energy transfer in the non- 
recoil approximation.) This latter formula can be used 
to derive a close expression for Ap x valid for arbitrarily- 
large, finite objects in the large impact parameter limit. 
In that case, only small w's contribute to F(u), due to 
the effective exponential cut-off imposed by the modi- 
fied Bessel functions Kq and K\. This means that only 
long wavelengths are relevant (to which the object ap- 
pears as small), so that it can be described by its static 
polarizability. Then, the u integral can be performed 
numerically to yield 



A Px = (5.55165 7 
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For comparison, the momentum transferred to a charge 
e at a distance b from the beam is Ap = — (2e 2 /6w)x. 

The large-6 limit given by Eq. JHJl is compared in Fig. 
121 with more detailed calculations that include higher- 
multipole moments, as described below. Also, the small 
particle limit of Eq. (J5J is discussed in Fig. [31 



B. Arbitrary size 

For larger particles or for close electron-particle en- 
counters, higher multipoles become relevant in the in- 
duced forces Then, it is convenient to express the 
evanescent field of the electron in terms of multipoles 
centered at the particle, so that the external electric 
and magnetic fields admit the following decomposition 
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where L — (l,m), k = ui/c, j'i(fcr) = 1 1 ji{kr)Y L {v), L = 
— \hv x V is the orbital angular momentum operator, and 
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FIG. 2: (color online). Momentum transfer to small spherical 
particles by a passing 200-keV electron as a function of the 
distance from the trajectory to the center of the spheres b. 
The momentum transfer has been scaled using the velocity 
v — 0.7c, the sphere radius a, and the impact parameter 
b. The perpendicular component of the momentum transfer 
with respect to the trajectory Ap x (solid curves) has been 
represented for spheres of radius a — 10 nm, 50 nm, 200 nm, 
and 500 nm (notice the rapid increase in Ap x near b = a). The 
parallel component Ap z (dashed curves) is only shown for a — 
200 nm and 500 nm. Dielectric alumina spheres and metallic 
silver spheres are considered (left and right plot, respectively), 
respectively. The large b limit for perpendicular momentum 
transfer [Eq. JSJ] is shown by horizontal dotted lines. 



with 



(for v = E, M) are multipole coefficients given by 
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In particular, L = (I, m) is conserved for spherical par- 
ticles and one has a linear dependence ijf£ in = t\ ip^** , 
where t" are scattering matrices that are given by ana- 
lytical expressions in the case of ho mog eneous particles 
of dielectric function e and radius a |l5|: 

-ji(Po)Pij[(Pi) + Pofi(Po)ji(pi) 
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where po = ka, p\ — po\fe with Im{pi} > 0, and the 
prime denotes differentiation with respect to p and p\. 

At this point, it is convenient to write the operators L 
and (l/fc)V in matrix form. One finds 
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and the x and y components of N are obtained from (jHJ) 
by rotating the reference frame using rotation matrices 
for spherical harmonics |19|. Exactly the same matri- 
ces as above apply to L and (l/fc)V acting on Hankcl 

functions . Furthermore, these matrices satisfy the 
properties L + = L and N + = — N. 

Now, the electric field admits an expansion of the form 
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C± = y/(l±m+l)(lTm). 



Here, s m = 1 if m > 0, s m = (— l) m if m < 0, and Cm is 
the Gegenbauer polynomial |l8j . The impact parameter 
b is defined in Fig. ^ 

The induced field around the particle is given by simi- 
lar expressions obtained by substituting i\) v £ oxt by new co- 
efficients ip u £ md , and ji by the Hankel function hi [l9j . 



where the coefficients 
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are obtained from the above expressions. Similar formu- 
las are obtained for H cxt and for the induced fields E md 
and H md in terms of multipole coefficients. Finally, we 
insert them into Eq. Q and perform the integral over a 
sphere in the s — > 00 limit. Then, the first two terms in- 
side the integrand give a vanishing contribution because 
the induced far-field is transverse. The remaining part of 
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the integral can be recast, noticing that only real terms 
must be retained, 
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andn(ft) = y/toFfi [(x+iy) Y 1 _ 1 /V2- (x-iy) Y n /V2+ 
z Yio] is the radial vector as a function of the polar di- 
rection Q. 



III. RESULTS AND DISCUSSION 

Fig. |21 shows the dependence of the momentum trans- 
fer on electron impact parameter b for alumina and sil- 
ver spheres of different radius, as calculated from Eqs. 
and © . Measured optical data have been used for the di- 
electric function of these materials 7] . One can observe a 
nearly exponential decay of the momentum transfer with 
b. Besides, the momentum transferred along the direction 
of the electron velocity vector (Ap z , dashed curves) is 
generally smaller than the remaining perpendicular com- 
ponent (Ap Xl solid curves), which finds an explanation 
in the fact that the contribution of these components to 
the energy loss Tito is vAp z + (Ap x ) 2 /m, where m is the 
electron mass: since mv 3> Ap, Ap x is allowed to take 
larger values than Ap z for each fixed lj. 

Notice also that Ap x converges quickly to the large b 
limit [Eq. H3 dotted curves] , producing a finite result un- 
der the scaling of Fig. |21 unlike Ap z , which goes faster to 
for large b. In this limit, the electron induces a dipole in 
the particle directed towards the electron, which results 
in an attractive force between these two similar to the 
image potential at surfaces |20j , leading to a momentum 
transfer Ap « Ap^x. For small metallic particles and 
closer encounters this picture is no longer valid and Ap x 
can actually reverse its sign and have a net repulsive be- 
haviour (e.g., in Fig. |21 for Ag particles of radius a = 10 
nm and also for the fullerenes of Fig.^J. 

A more detailed analysis of the magnitude of the mo- 
mentum transfer effect is given in Fig. EJ The momentum 
transfer is normalized to the particle mass M and the re- 
sult is the change in the particle velocity induced by the 
passage of the electron as a function of particle radius a. 
The trajectory of the 200-keV electron under considera- 
tion passes 10 nm away from the surface of the spherical 
alumina particles. The full-multipolc calculation [Eqs. 

and ©, solid curves] agrees well with the small par- 
ticle limit [Eqs. and dashed curves] when a is 
much smaller than b — a = 10 nm. Even though the 
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FIG. 3: (color online). Particle size dependence of the mo- 
mentum transfer normalized to the particle mass M under 
the same conditions as in Fig. |5] small particle limit (dashed 
curves) versus full multipole calculation (solid curves). The 
particle is made of AI2O3 (density p = 4.02 g/cm 3 ), the elec- 
tron energy is 200 keV, and the distance from the trajectory 
to the particle surface is 10 nm. Dotted curves show the mo- 
mentum transferred from light in an optical trap (see text for 
details). The inset depicts the change in the particle angular 
velocity as a result of the torque exerted by the electron. 



electron-particle interaction increases with the radius a, 
the actual change in the particle velocity shows a nearly 
exponential decay with increasing a. 

In a situation where the particle is trapped by lasers 
(e.g., in optical tweezers Q or in optical stretchers p|), 
one should compare the interaction with the electrons to 
the interaction with the laser light. To this end, we will 
consider a trapping cw-Ti:sapphire 100-mW laser emit- 
ting at a wavelength A = 785 nm and focused on a region 
of radius Rf = 10 /im. Furthermore, we will contem- 
plate the momentum transferred by the laser during the 
average time span At between two consecutive passing 
electrons in a transmission electron microscope operating 
at a current of 1 nA. The particle polarizability a is all 
that is needed to calculate light forces for the small radii 
under discussion (a A), according to Eq. Now, 
for real a this equation defines a conservative gradient 
force that responds to the potential — (a/2)|E| 2 , where 
E is the laser light field, whereas the imaginary part of a 
represents photon absor ptio n by the particle that trans- 
lates into light pressure 21]. These two components are 
represented separately in Fig. |21 after multiplication by 
At/M (dotted curves). The light pressure contribution 
is calculated for an incidence plane wave with the same 
photon flux as the laser at its focus. The gradient force 
component is obtained from the maximum force in the fo- 
cus region assuming a Gaussian profile for the laser field 
intensity (i.e., |E| 2 oc exp[-R 2 /(R f / In 2) 2 ]). Finally, it 
is convenient to define the polarizability from its relation 
to the scattering matrix, which upon inspection permits 
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writing a = (3/2fc 3 )tf. Unlike the well-known expression 
[l2| a = a 3 (e — l)/(e + 2), the former relation predicts 
a non-zero value for Im{a} even for particles with real 
e (like our alumina spheres), arising as a pure retarda- 
tion correction associated to radiation scattering (this is 
actually the origin of the light pressure component of 
Fig. I3J) . (Incidentally, gravity would produce a velocity 
change gAt = 1.56 nm/s, which is well compensated for 
in currently available optical trapping systems.) 

An important conclusion that can be extracted from 
Fig. Ellis that the crossover of trapping light into the main 
source of momentum occurs for particles of 20 nm in di- 
ameter when the electrons pass at a distance of 10 nm 
from the particles surface, thus allowing one to perform 
energy loss analysis of the transmitted electrons with sig- 
nificant statistics. Therefore, transmission electron mi- 
croscopy can be combined with in- vacuo optical trapping 
to study particles of sizes above some tens nm. 

While the transfer of momentum by the trapping light 
occurs in a continuous smooth fashion, the electrons de- 
posit all of the momentum during a small time interval 
~ a/v (<C At = 0.16 ns for 1 nA electron current). How- 
ever, the change in particle velocity per electron (vertical 
scale in Fig. |3|) produces a minute particle displacement 
during At (smaller than 1.6 x 10~ 9 nm -C a), and there- 
fore, the effect of the passing electrons is experienced by 
the particle as a nearly continuous source of momentum 
that is describable by an average force Ap/Ai. Actually, 
Fig. 01 suggests that using more intense electron beams 
(with even smaller impact parameters) acting during pe- 
riods of the order of one second will still not produce 
ejection of the particles from their trapping locations. 

It should be stressed that the momentum transfers that 
we have calculated using classical electromagnetic the- 
ory must be understood as the average value over many 
incoming electrons, since the actual strength of the in- 
teraction is not large enough as to guarantee that many 
photons are exchanged between each electron and a given 
particle. Like in aloof EELS experiments most elec- 
trons will not interact with the particles at all, so that the 
present results must be understood under the perfectly 
valid perspective of an statistical average performed over 
many beam electrons. The quadratic deviation from 
these average forces can play also a role (similar to strag- 
gling in stopping power theory), but this subject is left 
for future consideration. 

We have also studied momentum transfer to Ceo clus- 
ters (Fig. 01). The scattering matrices t\ have been ob- 
tained within the discrete-dipole approximation |22l |23| , 
where each carbon atom is described by an induced dipole 
whose polarizability is fitted to reproduce correctly the 
measured optical response of graphite 0. Further de- 
tails concerning the procedure followed to obtain t" will 
be given elsewhere j24j. At relatively small interaction 
distances b, the z component of the momentum is larger 
than the x component and the latter is negative. These 
are effects that can be hardly found in the above exam- 
ples and that originate in high-order multipoles (actually, 
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FIG. 4: (color online). Momentum transferred from a 5-keV 
electron to a Ceo cluster as a function of impact parameter b. 
The momentum is normalized to the cluster mass M. 

I < 5 are needed for convergence within the range of b 
under consideration). Even at a distance of 9 nm (no- 
tice that Ceo has a diameter of only 0.7 nm) the change 
in velocity produced by the passing electron can be sub- 
stantial. Therefore, the interaction of fast electrons with 
small clusters can produced dramatic effects if these are 
not mightily bound by a mechanism stronger than optical 
trapping. 

Finally, the passing electron can induce a torque on 
the particle that changes its angular momentum (AL y ) 
and makes it rotate. This is the effect discussed in the 
inset of Fig. |3 which shows the change in angular veloc- 
ity per electron, AJ7 = AL y /I, where / = (2/3)a 2 M is 
the moment of inertia of the alumina sphere. Like the 
electromagnetic force above, the torque is obtained from 
the integral of Maxwell's stress tensor ^3] i an d the de- 
tails follow a similar derivation as the one presented in 
Sec. [n] Averaging over the electrons of a 1 nA electron 
beam passing at 10 nm from the surface of an alumina 
sphere of radius a = 20 nm, one finds an angular accel- 
eration of 39 MHz/s. Under these conditions, the linear 
momentum transferred by the electrons can be absorbed 
by the trapping light, as discussed above. However, the 
angular momentum is not absorbed, and the particle will 
spin with increasing angular velocity until either the cen- 
trifugal force breaks it apart or radiation emission at the 
rotation frequency (vacuum friction) compensates for the 
electron-induced torque. 

In conclusion, we have shown that fast electrons follow- 
ing aloof trajectories (i.e., without direct overlap with the 
sample) in a transmission electron microscope can exert 
time-averaged forces on small particles of similar mag- 
nitude as those forces associated to trapping in optical 
tweezers and stretchers, and therefore, this effect can be 
used for analytical studies of mechanical properties of 
such particles, while electron energy loss spectra can be 
actually taken without causing ejection of the particles 
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from their trapping positions. 
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